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Abstract. One of the most surprising discoveries in quantum chaos was that nodal 
domains of eigenfunctions of quantum-chaotic billiards and maps in the semi-classical 
CN ' limit display critical percolation. Here we extend these studies to the level sets of the 

OQ , adjacency eigenvectors of d-regular graphs. Numerical computations show that the 

C*~) ■ statistics of the largest level sets (the maximal connected components of the graph for 



which the eigenvector exceeds a prescribed value) depend critically on the level. The 
critical level is a function of the eigenvalue and the degree d. To explain the observed 
behavior we study a random Gaussian waves ensemble over the d-regular tree. For this 
model, we prove the existence of a critical threshold. Using the local tree property of 
5_i ■ ci-regular graphs, and assuming the (local) applicability of the random waves model 



we can compute the critical percolation level and reproduce the numerical simulations. 
These results support the random- waves model for random regular graphs, suggested in 
PQ and provides an extension to Bogomolny's percolation model |2] for two-dimensional 
chaotic billiards. 



1. Introduction 

The statistics of the adjacency spectrum of random <i-regular graphs (in the limit of large 
vertex number) displays the generic attributes associated with the spectra of quantum 
Hamiltonians (in the limit h — > 0) whose classical dynamics is chaotic. This observation 
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which was first based on numerical simulations [3] , was further substantiated in [U [5] 
where a link was made between the distribution of cycle counts in ci-regular graphs, 
spectral statistics and Random Matrix Theory. In a way, the association of large d- 
regular graphs with chaotic dynamics is natural: The time evolution of random walks 
on a typical regular graph is known to mix exponentially fast [61 [7J. At the same time, 
the adjacency operator of a random regular graph is a symmetric matrix with identically 
(though not independently) distributed random variables (see section 12^21) . Thus, it does 
not come as a surprise that in the limit of large graphs with a fixed degree, the spectral 
statistics of a typical regular graph should follow some of the universality classes related 
to systems with chaotic (or mixing) dynamics. 

Recently, the study of level sets of eigenf unctions for two dimensional billiards 
and in particular the zero level sets (which are the nodal domains) have gained a 
considerable attention as possible indicators for the dynamics of the underlying classical 
system [H [9j [21 [101 HU [12]. The distribution of the (normalized) number of nodal 
domains in two dimensional billiards was examined in [9\. For chaotic billiards, this 
distribution was observed to converge in the semi classical limit into a universal measure, 
independent of the idiosyncratic dynamical characteristics of the billiard. These findings 
found an intriguing explanation in [2], where it was suggested that the distribution of 
nodal domains can be approximated by a non-correlated percolation process. Further 
connections with percolation theory were discovered in (TTJ [12] who showed that the 
boundary of the percolating domain reveal the SLEq statistics. These observations 
motivated the research reported in the present paper, where we examine the distribution 
of the adjacency eigenvectors for typical regular graphs. In particular we investigate the 
morphology of the associated level sets: Given a graph G and a real function f(G) which 
is defined on the vertices of G, the a-level sets of / in G are the maximal connected 
components of G, on which / exceeds the value a. 

Level sets of a special interest are the zero level sets - the nodal domains. 
In a previous paper [I], the dependence of the expected number of nodal domains 
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on the spectral parameter was studied. The numerically observed patterns were 
accurately reproduced by assuming a random waves model where the distribution of 
the eigenvectors converges to that of a Gaussian random field on the regular tree T^. 

Here we study further the morphology of level sets, and in particular investigate 
the percolation transition observed numerically at a critical level a c which depends on 
the eigenvalue under consideration. Invoking again the random wave hypothesis we 
are able to reproduce the numerical simulations, and thus bring further supporting 
evidence for the applicability of the random waves model. In this way we provide 
the discrete analogue to Bogomolny's percolation hypothesis [2], and shed more light 
on the surprising connection between percolation and spectral theory. While for two 
dimensional billiards, if a critical level set exist it must be the nodal set by duality 
arguments (e.g. p3]), the fact that in the present case the percolation threshold typically 
occurs at a c ^ should not come as a surprise. This is consistent with classical results 
about probabilistic percolation on graphs [14]. 

The paper is organized in the following way. After some necessary preliminaries and 
a review of some pertinent results, we describe the results of numerical simulations where 
the distribution of the largest a-level set for different eigenvectors of random regular 
graphs were studied. The numerically deduced critical level a c (X,d) as a function of 
the eigenvalue and the degree of the graph summarizes these computation. The next 
chapter is dedicated to a systematic construction of a Gaussian random waves model 
on the d- regular tree [15]. For this process, we are able to prove that a percolation 
transition occurs. We then show that the critical level sets computed for T^ reproduce 
accurately the numerical simulations for d regular graphs, which is reasonable in light of 
the local tree property. Yet, a rigorous proof of this observation is lacking. The paper is 
concluded by a comparison between the observed critical behavior and the percolation 
model for two dimensional billiards. 
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2. Preliminaries 

2.1. Elementary definitions 

A graph G is a discrete set of vertices, connected by edges. We denote the size of a graph 
by \G\ = n, where by an abuse of notation, we use the symbol G to denote both the 
graph and its set of vertices. We will consider only simple graphs, i.e. graphs containing 
no loops or multiple edges. A graph G is d-regular if for every vertex v G G, the degree 
of v (or the number of edges connected to v) is exactly d. For vertices v, v' G G, we 
define the graph distance \v — v'\ as the length of the shortest walk in G from v to v'. 

A graph is completely specified by its adjacency operator A, where A^ = 1 if 
Vi, Vj G G are adjacent vertices, or zero otherwise. As the adjacency operator is real and 
symmetric, it has n real eigenvalues. We will denote the spectrum and eigenvectors of 
the adjacency operator by 

Ai<A 2 ... <A n (A-X J )f j = 

The eigenvectors are normalized through the paper according to the convention 

£/» = » 

veG 
so that Var(/(v)) does not vanish as n — > oo. 

2.2. G(n,d) - the ensemble of random regular graphs 

For a given value of n, d G N, the ensemble G(n, d) consists of all d regular graphs on n 
vertices, equipped with the uniform measure. 

The geometrical and spectral properties of G(n, d) have been extensively studied for 
more than 30 years, and are successfully applied in various fields such as combinatorics, 
information theory, pseudo-randomness and more (see [16] for a review) . In the following 
we would like to investigate the properties of the eigenvectors of a typical (n, d) graph, 
for a fixed d > 3 in the limit n — > oo. Here and in the following, by stating that a typical 
(n, d) graph has a property X, we mean that the probability that a graph G G G(n, d) 
has the property X converges to one as n — > oo. 
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In [T7] the distribution of short cycles in G(n, d) is calculated. Denoting by C& 

the number of independent A;-cycles in a graph, it is shown that the distribution of 

{Ck}, 3 < kle\og d _ 1 n/2 converges as n — > oo into independent Poisson random 

variables with an expectation value 

id - l) k 

nc k ) = { —^- . (2.i) 

Since the expected number of fc-cycles does not increase with the size of the graph, 
we find that a ball of radius c log n around a random vertex has a tree structure with 
probability 1 — n c_1 , for every c < 1/2. As a result, the local structure of an (n, d) graph 
near most of its vertices is identical to that of Td - the <i-regular (infinite) tree. 
The diameter of an (n, d) graph, i.e. the maximal distance between vertices in G, is 
given by [T8] 

diam(G) = log(nlogn) + 0(1) (2.2) 

(here and in what follows the logarithm base is (d — 1)). This result shows that the 
typical distance between vertices along the boundary of the 'local tree' is of the same 
magnitude as the distance between two arbitrary vertices in G. 

The local resemblance between (n, d) graphs and Td is reflected in the spectral 
density of their adjacency operator as well. Both the spectrum of the tree [19] and the 
limiting spectral density of G(n, d) [20] are supported on the interval 



a(T d ) = \-2Vd-l, 2Vd - 1] (2.3) 

with a spectral density, given by 



For every connected regular graph which is not bipartite, the unique stationary 
distribution for random walks on the graph is uniform over the vertices. The rate of 
convergence into the stationary distribution is dictated by P^J, - the probability that 
a walk of length k which begins at the vertex v will terminate in v'. In [7] it is shown 
that for a typical (n, d) graph 



lim sup 

k— ¥oo 



p{k) _ 1 
v,V n 



l/k 

= e~ 7 (2.5) 
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where the Lyapunov exponent 



1 2y/^T „( 1 
7=1 : YU 



d \ log n 

is strictly positive for d > 3. As a result, random walks on a typical (n,d) graph are 



exponentially mixing with a Lyapunov exponent 7^1 — 2y/d — 1/d. This observation 
justifies attributing the title "chaotic" to (n, d) graphs, as was done in the introduction 
section. For further detail on expanding (equivalently, mixing) graphs and the relations 
between their spectral, geometrical and dynamical properties we refer the reader to [16J. 

2.3. The random waves model for the adjacency eigenvectors of G G G(n,d) 

The main tool in the present study is the random wave model [H [15] . It is based on the 
observation that for the d regular tree Tj, the distribution of a typical eigenvector can 
be approximated by a real Gaussian process GSd{X). The process associates random 
functions i/) u : Ta — > 1R to the regular tree, so that for every subset of vertices {vj} C Td, 
the distribution of any linear combination of ipio{ v j) is Gaussian. 

A Gaussian process is uniquely characterized by its mean and covariance operator. 
Introducing the Chebyshev Polynomials of the second kind 

= an((t + l)coB-(»)) 

sin (cos -1 (a;)) 



and the polynomials 119, 21] 

It was shown in [TJ [15] that for every A G o-{T d ) (12.31) . the random Gaussian process 
QSd{X) which is characterized by the covariance 

Cov(Mv),Mv')) = W (\v - v'\) (2.8) 

Has the following properties: 

(i) (Ar d — XI)ip u — for almost every ip^ G QSd(X)- 
(ii) QSd(X) is invariant with respect to the symmetries of T^. 
(iii) The process is normalized, so that V&r(ip u (v)) = 1. 
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The random wave model is based on the conjecture that the distribution of a typical 
adjacency eigenvector of a graph G G G(n, d) graph, with an eigenvalue A G cr{Td) 
converges locally to that of QS^iX). This hypothesis, was supported by various numerical 
tests, and found a partial formal justification in [15]. The random waves hypothesis will 
be the basis for the analysis of the morphology of level sets in (n, d) graphs, which will 
be carried out in this paper. It is the analogue of Berry's model for the distribution of 
eigenfunctions for chaotic billiards |22j . 

3. Level sets percolation on (n, d) graphs 

In this section we shall present the numerical evidence which led us to propose that 
level sets undergo a percolation transition in the limit n — Y oo. 

For a graph G, a real function f(G) and a given a G ffi, we denote by 

G a {f) = {veGJ{v)>a} (3.1) 

the induced graph, which is obtained by deleting all the vertices for which f(v) is below 
the threshold a; The a-level sets of / in G are the connected components of G a (f). 

Since there is no known analytical expression for the distribution of eigenvectors in 
an (n, d) graph, we cannot offer an expression for the limiting distribution of the a-level 
sets for this ensemble. However, motivated by the resemblance between the spectral and 
eigenfunctions statistics for G(n, d) and chaotic billiards, we have looked for a numerical 
evidence to a phase transition in (n, d) level sets. 

For a given graph G G G(n,d), a G R and an eigenvector (A — XI) f = 0, we define 
Ga to be the largest component of G a (f) (13. II) and evaluate the ratio \G™ \/\G a \. 
We have generated (following [23]) and diagonalized (using MatLab) random regular 
graphs on up to 4000 vertices with degrees ranging from 3 to 15. For each eigenvector / 
of each graph we have measured the ratio \Ga |/|Cr a | while varying a from min(/) to 
max(/). The ratios \G a \/\Ga | are plotted in figure [T^s a function of a for a 3-regular 
graph on 4000 vertices. The different lines correspond to different values of A are given 
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in the inset. 
A sharp transition in the normalized size of the level sets is evident: for every value 




Figure 1. The ratio between the magnitude of the largest level set to the size of the 
induced graph G a for a single realization of a (4000, 3) graph. Each curve corresponds 
to one eigenvector, whose level sets are interrogated by increasing a. 



of A there is a narrow window in the vicinity of some a c (X,d), so that for a > a c , 
the ratio \Ga \/\G a \ is close to zero, while for a < a c , \Ga \/\G a \ is of order one. 
The described phenomenon was observed for all the tested values of d and for all the 
examined eigenvectors. 

Moreover, repeating the experiment, while varying the size of the graph, we have 
observed that the value of a c (X, d) does not vary with n, while the transition becomes 



sharper as n increase. As an example, in figure |2] we plot the variation in \G. 



(max) I 



i a 



for the eigenvectors which correspond to A ~ 0, for different realizations of 3- regular 
graphs of various sizes. 

This extensive corpus of numerical data provides strong evidence supporting the 
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Figure 2. A comparison of the curve |Gq |/|G Q |, for several 3-regular graphs of 
varying size, the curves correspond to n = 100, 250, 1000, 4000 where in each graph we 
consider the eigenvector with the closest eigenvalue to zero. 



existence of a phase transition for the level sets of (n, d) eigenvectors. Namely, for every 
d > 3 and A G o~(T ( i), there may exist an a c (X, d) G M, so that the level sets of a typical 
(n, d) eigenvector which corresponds to the eigenvalue A are all microscopic for a > a c , 
while for a < a c a macroscopic component is expected to appear. 

Note that the suggested transition differs from the percolation hypothesis for chaotic 
billiards in two main aspects. First, while for billiards the transition is expected to follow 
the characteristics of non-correlated percolation, for (n, d) graphs we expect correlations 
to be relevant (as will be discussed in section [5]). Second, Unlike the percolation model 
for billiards, the critical threshold for regular graphs depends on the corresponding 
eigenvalue. 

In order to estimate the dependence of a c (X,d) on its arguments, we have chosen 
(somewhat arbitrarily) to identify a c (A, d) with the steepest point of the curve 
\Ga \/\G a \, obtained for a graph of size n = 4000. In figure [3] we present our numerical 
estimate of the critical curves. The results suggest that for a given value of d, the critical 



Percolating level sets of the adjacency eigenvectors of d-regular graphs 



10 




-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 



2 v / d Zr T 



Figure 3. A numerical estimation to a c (A, d) for d = 3,5, 6, 12. 

threshold increases monotonically with A for all d > 3. For the lowest degree, d = 3, 
the critical curve a c (A, d — 3) shows a minimum in the vicinity of A ~ —0.52. 
In the next sections we shall show that the numerical resulted summarized above can be 
reproduced theoretically by considering first the level sets in the random process QSd{\) 
on Trf, and then assuming the validity of the random waves conjecture for the adjacency 
eigenvectors in G(n, d). 

4. The distribution of level sets in QSd(X) 



In this section we shall study the distribution of the level sets in GS and prove that 
they undergo a percolation transition, for which the critical threshold a c (X, d) can be 
computed. 

As the process QSd{X) is Gaussian and characterized by the covariance operator 
(12. 8p . it is possible to rigorously analyze its level sets statistics. Setting 

T u (a) = {veT d} if> u {v)>a} (4.1) 
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to denote the induced a level sets tree for a given ip^ G QSd(X), we show that: 

Theorem 4.1. VA G o~(Td), there exists an a c G R so that for almost every realization 
ipu £ GSd(X), T u (a) has an infinite component for a < a c , but only finite components 
for a > a c . 

As the proof of the theorem is rather technical, we refer the interested reader to 
[T5] where a complete and detailed proof of the theorem can be found. Here, we shall 
provide the main line of the proof, skipping much of the technical aspects. 

It is important to note that for a given a G R, the a-level sets of the process QSd(X) 
is a homogeneous vertex process on Td- That is, the probability measure of the process 
is invariant with respect to the symmetries of Tj. 

For any homogenous process, the probability that v, v' G Td belongs to the same 
connected component depends only on the distance between the vertices and will be 
denoted by P\ v - V '\- Denoting the sphere of radius k around v G T d by 

A k (v) = {v' ET d ,\v-v'\ = k} (4.2) 

we find that the probability that v is connected to its fc-sphere is at most |A&(i;)|P&. 
Since |A fe | = d(d — l) fc , we obtain that if 

limsu P P, 1/fc <-^- (4.3) 

fc— >oo a l 

then the probability that the connected component of v exceeds the radius k decays 
exponentially with k, so that the probability to find an infinite component is zero and 
the process is subcritical. 

We should note that the opposite statement is not necessarily correct, i.e. there are 
percolation processes on Td for which lim inf ^oo P k > l/(d — 1), but do not contain 
infinite components, due to long range correlations. 

The influence of the correlations which are induced by the process can be formally 
evaluated according to the next classification of random processes on trees, introduced 



in 
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v eC ffl (v 0AlA2 ) 



Figure 4. For a quasi-bernoulli process, the probability to find vo € C%(v±) 
(i.e. that the continuous purple line is entirely occupied in T u ) conditioned that 
vq 6 C"{v%) (dotted green line is occupied) is uniformly bounded by the probability 
that vo € C"(vi), conditioned that vo 6 C"(v A v\ A Vo) (dashed blue line is 
occupied). Therefore the effect of values which are obtained by the process along 
the path i>oaia2 — > V2 has only a limited effect on the path «oaia2 — > v%. 

Definition 4.2. A random percolation process on a tree graph T, associating for all 
io G Q an induced subgraph T w C T, is a quasi Bernoulli process, if 3M < oo, such that 

Wv ,v 1 ,v 2 e T: 

F ( Vl e C r ^v )\v 2 e C^(v )) ^ 
P(«i G Cr w (^o)koAiA2 G C Tu (v )) - { ' ] 

where v oaia2 is the intersection of the simple paths in T between the three vertices (see 
figure HJ) and Cr„(vo) is the connected component of v in T u . 



According to definition 14.21 in a quasi Bernoulli process, the distribution of the 
process along any simple path is only moderately affected by the realization of the 
process outside the path. As an example, it can be verified that a Markovian process is 
quasi Bernoulli. 

The quasi Bernoulli classification provides a simple criterion for the relevance of the 
long range correlations of the process on its macroscopic properties, as suggested by the 
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next lemma [24]J: 

Lemma 4.3. (Lyons) Let {Q, P} be an invariant quasi Bernoulli process on Td, which 
associates Vw £ Q an induced graph T w C T^. If\§ 

lim P l k ' k < -J- (4.5) 

fc->oo a — 1 

then, all the connected components of T u are almost surely finite. If 

Jim Pl ,k > -±- (4.6) 

k—>oo a — 1 

T w will almost surely have a component of an infinite cardinality. □ 

An equivalent phrasing of the lemma is the following: for a quasi Bernoulli process 
on Td, all the components are almost surely finite if their expected cardinality is finite, 
while if the expected cardinality diverges, infinite components will almost surely exist. 
As the a- level sets are monotonically decreasing in a, we find that in order to prove 
theorem 14.11 it is enough to show that 

• For every A G cr(Td) and a6l, the a-level sets of QSd(X) are quasi Bernoulli 

• For small enough values of a there almost surely exist infinite level sets, while for 
large enough values the a-level sets are all finite with probability one. 

The first condition guaranties that for every A G o~(Td) and a G M, the probability to 
find an infinite component is either zero or one (according to lemma H~B1 . Assuming the 
level sets are quasi Bernoulli, and since lim^oo Pj (ex) is monotone in a, the second 
condition ensures that for every A G R a supercritical and a subcritical phases exist, 
where the transition between the two occurs at a c (X,d) which is given by the implicit 
expression 

lim Pl /k (a c ) = J- (4.7) 

k— >oo a — 1 

| The theorem as it appears in [24j characterizes general quasi-Bernoulli process on an arbitrary 

tree graph. For the sake of clarity, we consider here only the restriction of the theorem to invariant 

percolation processes on Td- 

§ Note that for a quasi Bernoulli process, the limit lim| 11 /_ u |^ 00 (P(i/ € Ct u (v))) '' v ~ v ' does exists. 

This can be verified by restricting the condition (|4.4|l to cases in which vo is on the simple path between 

«i to i> 2 , i.e. v = w aia2- 
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To verify the existence of a subcritical regime, we note that if a simple path U of 
length k is contained in an a-level set, then necessarily \/vj G U, ip u (vj) > a. Therefore, 
setting ^ui(U) = Y1j=i' 1 Poj( v j)j we obtain that 

P k {a) = F(y Vj 6 U, ^ u {vj) > a) (4.8) 

<P(* W (E7) > ka) 
Note that ^^(U) is a Gaussian random variable, with mean zero and variance 



Var(^([/)) = E 



5^"K)VU%) 



(4.9) 



= fc (V A) (0) + 2^^0 {A) (j)) < fc$W 

where $W = 0^(0) + 2J2f =1 \<f>M(j)\ ( 12~7|) . 

Since (f>^ x '(j) is exponentially decreasing in j (12.71) . we find that $^ A ^ < oo. Therefore 

Psw< v»r exp (-^) <e "" (4io) 

where /3 = (2$W) -1 . 

As |Afc| = d(<i— 1) x , we obtain from (14. 31) and (14.101) that the a- level sets are subcritical 

for a > ^/ln(d— l)//3. 

The existence of infinite a-level sets for small enough a is a straight-forward 
consequence of [25J, where it is shown that every homogenous vertex percolation process 
on Tj, for which the survival probability exceeds d/2(d — 1) is supercritical. As ipu(v) 
distributes as a normal variable with mean zero and variance one, the vertex survival 
probability for T a {%l) u ) is given by 



CO 



Po(o) = -i/ e- 2 / 2 . (4.11) 

V 2,71 Ja 

Since d/2(d — 1) < 3/4 for every d > 3 and as according to (14. lip Po(a) > 3/4 for 
a < —0.68, we obtain that for every d > 3 and A e cr(Td), the a-level sets of £?<S<f(A) are 
supercritical for a < —0.68. 

In order to prove that the a-level sets of GSd(\) are quasi Bernoulli, we have to 
show that the long range correlations do not dominate the structure of the random tree 
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T u (a) (14. ip . A major step toward this goal is the next theorem, which identify the 
following Markov property of the underlying process QSd(\)'- 




Figure 5. A set of vertices which corresponds to the requirements of theorem 
14.41 t>o,vi an d i>2 are denoted explicitly; {vj}j>3 are marked by full red circles; 
{V 6 Afc, \v' — V\ | < \v' — V2 1} are marked by hollow black circles fsee l4.12]) for the case 
fe = 3. 



Theorem 4.4. Let {vj}j =0 C Td (for k > 3), so that \v\ — v%\ = 1, and V\ is on the 
simple path between vq to {fj}j =2 (see figure^). Then VA G o-{Td), the distribution of 
4>u){vq) conditioned on {ip u (vj)}j =1 is independent of {ipuj(vj)}^ =3 

The proof of theorem l4.4l relies on the observation that for any adjacency eigenvector 
ip : Td — > R, the value of ip(v\), ^(^2) determines the value of the sum 

£ W) (4-12) 

v'EAk(vi),\v'—vi\<\v'—V2\ 

(see figure [5]). Since the distribution of QSd(X) is invariant with respect to the 
symmetries of Td, we obtain from the last observation that given {if} u (vj)}j =1 , the 
expected value of ^{vq) is determined by ip u (vi), ^(^2)- 

Since non-correlated components of a multi-normal random vector are also independent 
(e.g. [26J), the theorem follows. □ 
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The described dominance of short range correlations in the process QS^iX), is 
utilized in [15], to bound rigorously the effect of long range correlation in T w (a), resulting 
in the identification of T w (a) as a quasi Bernoulli process and in the establishment of 
theorem 14.11 

The last step in the analysis is obtained by solving (14. 7p for various values of d 
and A G er(Td). The results are shown in figure [61 superimposed on the numerical 
data obtained for G(n,d) graphs and discussed in the previous section (figure [3]). The 
agreement between the G(n, d) data and the critical threshold levels computed for the 
QSd{X) process on Td is perfect (including in particular the non monotonic behavior 
of a c (\,d = 3) which is reproduced as well). It strongly supports the random waves 
conjecture for G(n,d). As was already stated above, the random wave conjecture is 
valid locally. In the present context, however it applies globally. The proof of the 
random wave conjecture is still lacking and is a challenge for experts in probabilistic 
graph theory. 

5. Discussion 

The identification of a phase transition for the level sets of QSd{X) (and the observed 
transition for G(n, d) eigenvectors) may be considered as an extension of the percolation 
hypothesis for random waves on IR 2 and for chaotic billiards [2]. We would like to 
conclude this paper by a short comparison between the two models. 

We should note that the formal derivation for both of the model is based on the 
statistical properties of the corresponding random waves model, while the heuristic 
justification of the percolation hypothesis for two-dimensional random waves is based 
on arguments which are not clearly justified [271 I2H1 [2U E2], theorem 14.11 provides a 
rigorous proof to the existence of a critical level-set for QSd(X)- 

In both of the cases, the applicability of the suggested model to the corresponding chaotic 
system (a billiard or a mixing regular graph) is a consequence of the corresponding 
random waves model and not of an independent derivation. However the consistency of 
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Figure 6. A comparison between the value of a c (X,d) for the Gaussian process 
GSd(^), given by equation I4.7I (lines) to the numerical estimation of a c (A, d) for (n, d) 
graphs (markers). 



the models with various numerical tests, such as [9j |2j [30l ED] (for billiards) and [T] (for 
(n, d) graphs) provide a firm support both for the percolation and the random-waves 
models. 

The arguments behind the two percolation models differ significantly. The 
percolation hypothesis for two dimensional chaotic billiards is justified in [2j [10] by 
dimensional arguments which prevents its generalization for generic chaotic systems, 
or even billiards of higher dimensions. In particular it is based on the topological 
identification of a random wave with the square lattice [2J [12j [28] and the particular 
critical threshold p c = 1/2 for self-dual bond percolation processes on the square lattice. 
In addition, the neglect of correlation is based on a careful application of the Harris 
criterion [31J, which is valid only for two-dimensional systems [T0[ 125]. 
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Unfortunately, the identification of a critical level set for the process QSd{\) 
demonstrates the same weakness, as it utilizes repeatedly the tree structure of T^. 
Nevertheless, the identification of a critical level sets for the two seemingly non-related 
waves ensembles hints on a universal mechanism behind the phenomena. In addition 
the dependence of the critical point (J4.7P in the spectral parameter A, implies that 
the suggested transition for the eigenfunctions of a generic chaotic system (assuming it 
indeed exists) may exhibit a more complicated behavior then the simple model which 
is suggested in [2]. 
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